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SUMMARY We propose a new method to construct a four 
parameter family of quantum-mechanical point interactions in 
one dimension, which is known as all possible self-adjoint ex- 
tensions of the symmetric operator T = — A |"Cq° (R\{0}). It is 
achieved in the small distance limit of equally spaced three neigh- 
boring Dirac's <5 potentials. The strength for each <5 is appropri- 
ately renormalized according to the distance and it diverges, in 
general, in the small distance limit. The validity of our method 
is ensured by numerical calculations. In general cases except for 
usual <5, the wave function discontinuity appears around the in- 
teraction and one can observe such a tendency even at a finite 
distance level. 

key words: quantum mechanics, one dimension, generalized 
point interaction, wave function discontinuity, functional analy- 
sis 

1. Introduction 

Point interaction is a simple but useful object for ex- 
amining the effect of small impurities on low-energy 
dynamics in quantum mechanics. It serves as a mini- 
mal perturbation which gives light on the influence of 
more complicated perturbations on dynamical nature 
of the system. It has been first observed [1], [2] that 
wave chaos (chaotic nature in quantum spectrum) is 
induced by point perturbation in two-dimensional in- 
tegrable billiards and a general condition for its emer- 
gence has been clarified in [3]-[5]. In [6], the depen- 
dence of the "degree" of chaos on the number of point 
impurities has been examined in details. 

Historically, the first influential paper on point in- 
teractions was given by Kronig and Penney [7]. The 
Kronig-Penney model (potential consisting of a peri- 
odic array of S functions) has been widely considered 
as a standard reference model in solid state physics 
for more than six decades. In two and three dimen- 
sions, the point interaction is nothing but an extended 
object of the usual 6 potential. However the phe- 
nomena are much richer in one dimension. This is 
due to the fact that whereas the symmetric operator 
T = -A[C£°(R d \{0}) has deficiency indices (1,1) in 
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spatial dimension d — 2,3, it has deficiency indices 
(2,2) in dimension d = 1 [8]. As a result, the oper- 
ator T has in one dimension a four-parameter family 
of self-adjoint extensions [9]. The general connection 
condition is characterized by 

^ rt+o) ) ~ v \ <p(-o) ) ' [ ) 

where V takes a form 

V = e l8 U (2) 
with 9 6 R and 

W=(" J)eSL(2,R), (3) 

namely, a5 — (3"f = 1 [10]. Here we assume that the 
point interaction is placed at the origin on x-axis. The 
condition (1) covers the usual S potential, the connec- 
tion matrix of which is given by 

v=v«(«)=( j \ y (4) 

Here v is the potential strength. Another interesting ex- 
ample is a so-called (historically ill-called) 6' potential, 
which induces such a boundary condition that the wave 
function has continuous first derivative on the right and 
left, but it has a jump proportional to the first deriva- 
tive [11]; 

v=v.(«) 3 (i (5) 

where u is regarded as the strength of the interaction. 
It has been shown [12] that the connection condition 
(5) has little to do with the derivative of S potential. In 
this paper, we denote the contact interaction with the 
transfer matrix (5) by e{x). 

Except for S, the connection condition (1) appears 
at first sight to be unnatural for quantum mechanics, 
and some efforts have been made in order to seek the 
physical Hamiltonians which embody the condition (1). 
Seba has observed [13] that a certain two-parameter 
family is obtained by non-local pseudopotcntials with e 
potential. Chernoff and Hughes have shown [14] that a 
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three-parameter family disjoint from Scba's class corre- 
sponds to local pseudopotentials involving the so-called 
S 2 . Carreau has succeeded in realizing a general class 
in the small-size limit of non-symmetric operators [15]. 
A local realization has been given by Roman and Tar- 
rach by using a potential consisting of a step- like barrier 
sandwiched in two 6 potentials [16]. Their model re- 
quires, however, a subtle limiting procedure and indeed 
the height of the middle barrier as well as the strengths 
of the two side 5's are necessary to be renormalizcd 
with a doubly logarithmic dependence on the distance 
between the two <5's. In spite of considerable works, 
existing approximations do not seem to have much rel- 
evance to the experimentally realizable systems. Our 
main purpose is to rectify this situation. In [17], [18], 
we have constructed the e potential in the small dis- 
tance limit of three nearby S potentials and realized 
in terms of usual 5 and e potentials a three-parameter 
family of self-adjoint extensions under the assumption 
of time reversal symmetry. In this paper, we general- 
ize the previous formalism and give in a direct manner 
within the three <5's model a complete solution for the 
four-parameter extensions covering the cases without 
the time reversal symmetry 

In Sect. 2, we start by solving the one-dimensional 
Schrodingcr equation with a constant vector potential 
within the framework of the transfer matrix formal- 
ism. Considering the potential which consists of equally 
spaced three 5's with a constant vector potential be- 
tween the two side 5's, we show that all possible ex- 
tensions are attained in the short range limit of the 
assumed potential if the strengths of the <5's as well as 
that of vector potential are renormalized according to 
the distance in a simple manner. In Sect. 3, We give 
numerical examples which justify the method in the 
previous section. The current work is summarized in 
Sect.4. 

2. Approximation of Four-Parameter Family 
of One-Dimensional Point Interactions by 
Three Nearby Delta Interactions 

We first consider the Hamiltonian with a homogeneous 
(everywhere constant) vector potential in one dimen- 
sion. It is given by a usual minimal coupling to the 
magnetic field; 



H=(p-A) 2 



(6) 



where p = — id/dx is the momentum operator and con- 
stant A is the strength of the vector potential. The 
Schrodingcr equation for the Hamiltonian (6) reads 



d 2 ip(x) „ ,dip(x) , 9 , . 
^ v ' +2iA-^—^ + A 2 ip(x) 



fcV(x). (7) 



dx 2 ' dx 

Here k is the wave number of the particle. For later con- 
venience, we rewrite the equation (7) within the transfer 
matrix formalism. For this purpose, let us introduce 



a vector notation for the wave function and its space 
derivative; 



(8) 



In this notation, the equation (7) is rewritten by the 
first-order coupled equation 



¥'(ar) =H(A, fc)*(x) 



with 



H{A, k) 



2iA -k 2 + A 2 
1 



The solution of Eq.(9) is given by 
*(a;) = G(A, k;x- x )V{x Q ), 



(9) 



(10) 



(11) 



where G(A, k;x) is the exponential function of 
H(A,k)x; 



= p H(A,k)x 



e iAx |cos(fca;) ( J J ) 
sin(fca;) / iA -k 2 + A 2 



G(A,k;x) 



+ k \ 1 
Since 

TrTL{A, k) = 2iA, 
the matrix Q(A, k; x) satisfies 

detg(A,k;x) = e 2lAx , 



-iA 



■ (12) 



(13) 



(14) 



indicating that the matrix on the RHS of Eq.(12) has 
determinant one. 

In order to realize the general condition (1) in the 
small-size limit of a realistic finite-range potential, we 
consider three neighboring usual S's which are located 
with equal distance a [17]. In order to take into ac- 
count the time reversal symmetry breaking, we add a 
constant vector potential between the two side <5's; 



A a (x) = A9(a - x)9(x + a), 



(15) 



where 8 is a usual step function. The Hamiltonian is 
hence given by 



H a = (p-A a (x)) +v-6(x + a)+v 5(x) 
+v + 5(x — a), 



(16) 



where the strengths of the three delta's at x = ±a and 
are set to v± and vq, respectively. The assumed po- 
tential V a (x) is shown in Fig.l. Inserting Eq.(15) into 
Eq.(16), one obtains 

d 2 d 
dx z dx 



+ (w_ + iA)S(x + a) + v S(x) 
+ (v + — iA)8{x — a). 



(17) 
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V a (x) 



A a (x) = 



A a (x) = A 

v (a) 



—a 

"T 



A a (x) = 







X 



V-(a) v + (a) 



Fig.l Approximation of a general class of point interactions 
by three neighboring <5 potentials. The vector potential A a (x) 
has a step-like form: a constant A between two side 8's and 
otherwise zero. The strengths of the three 8's depend on 
the distance a in general. It will be shown that if i>-f (a) = 
-l/o+(a + l)/7, v (a) = 7/a 2 and v-(a) = -l/a+ (<5 + 
together with 6 = 2Aa constant, one attains in the small a 
limit the boundary condition (1) for 7 =|= 0. For 7 = 0, the con- 
dition is realized if v+ (a) = (a — I)/ (2a) , vo = 4/3/ (a + S + 2) 
and V- (a) = (<5 — l)/(2a), respectively. 



The appearance of the complex £'s at x — ±a are in- 
evitable for ensuring the Hamiltonian to be symmetric 
and they indeed describe a sudden change of the vec- 
tor potential at the location of the two side <5's. With 
the Hamiltonian (17), the connection of between 
x = — a — and x = a + is given by 

¥(a + 0) = V„(A,fc)¥(-a-0), (18) 
where 

V a (A,k) = V s (v+-iA)g(A,k;a)V s (vo) 

xg(A,k;a)V s (v_+iA). (19) 

Inserting Eqs.(4) and (12) into Eq.(19), we obtain 

V a (A,k) = e 2iAa U a (k), (20) 

where the (2, 1) component of U a {k) is given by 



\ii mi sin(2fca) , sin 2 (fca) 
[W (fc)J 2 i = ; h 



vo 



(21) 



121 k ' fc 2 

and the other elements are written as 
sin ( 2 K a } 

[U a (k)] n = cos(2ka) + '- v + [U a (k)} 21 v+, (22) 

[U a {k)] 22 = cos(2ka) + sm ^ a ) VQ + [ Ua {k)] 21 v., (23) 

[U a (k)} 12 = cos 2 (ka) ■ (v+ + v + V-) 
— sin 2 (ka) ■ (v + + V-) 
sin(2fca) 



+ 2k 

21 V + V 



{-2k 2 +v (v+ +v-)} 



in terms of [U a (A, k)] 21 , respectively. Note that the ef- 
fect of the magnetic field is absorbed in the phase factor 
in front of the matrix U a {k) in Eq.(20). Since 

detV„(A,jfe) = e 4tAa (25) 

from Eqs.(14) and (19), the real matrix U a {k) satisfies 

det U a (k) = 1, (26) 

namely, U a (k) G SL(2,H). 

The phase factor of the matrix (20) shows that in 
order to realize the general connection condition (1), 
one has to take a small-size limit under the condition 



2Aa = 6 



(27) 



with a fixed 9 e R. Accordingly, the strength of the 
vector potential increases as a — ► +0; 



A(a) 



2a 



(28) 



As shown just below, a singular behavior for small a is 
demanded for the strengths of the three <5's potentials 
in order to attain the condition (1). This enforces a 
careful analysis of each matrix element of lA a (k) in the 
small distance limit. Indeed, it requires the expansions 
of triangular functions appeared in Eqs.(21) - (24) up 
to the orders shown below; 



cos(2fca) = 1 + (9(a 2 ), 
sin(2fco) = 2ka - ^k 3 a 3 + (9(a 5 ), 



(29) 
(30) 

cos 2 (fca) = 1 - k 2 a 2 + 0(a 4 ), (31) 

sin 2 (fca) = k 2 a 2 - \k A a 4 + 0(a 6 ). (32) 

Inserting Eqs.(30) and (32) into Eq.(21), one observes 

[U a {k)] 21 =2a + v a 2 + 0(a 3 ). (33) 

This indicates that the strength vq of the middle S is 
required to diverge with the order of 1/a 2 in the small 
distance limit in order to make \U a (k)] 21 have a non- 
zero limit. Indeed, if u is changed according to the 
distance a as 

(34) 

with an arbitrary real constant 7 =|= 0, one obtains 

[U a {k)] 21 = 7 + 2a - ^k 2 a 2 + 0(a 3 ), (35) 
leading to 

a lrm o [W a (£0] 21 = 7 - (36) 
Using Eq.(35), one has the estimates 

[U a {k)] n = 1 + 1 + (7 + 2a + 0{a 2 )) v+ 



7 

v o( a ) = -o 
a 1 



Pa{k)] 



22 



(24) 



+0(a), 

^ + l + (7 + 2a+0(a 2 ))«- 
+0(a) 



(37) 



(38) 
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for small a. The RHS of Eqs.(37) and (38) diverge in 
the small a limit if one keeps v + and V- constant. The 
strengths v + , w_ of the side <5's are demanded to diverge 
with the order of 1/a to make [U a (k)] 11 and [U a (k)] 22 
converge in the small a limit; If 



v+(a) 
v-(a) 



1 a + 1 

- + , 

a 7 

1 (5+1 

a 7 



(39) 
(40) 



with real constants a and 5, the diagonal elements con- 
verge into 



lira . [W (fc)]n 

a »+0 

lirn [U a (k)] 22 

a >+0 



(41) 
(42) 



respectively. Note that the divergent term in v+ and 
u_ is dependent only on the distance a between neigh- 
boring (5's and the regular part determines the bound- 
ary condition around the generalized point interaction. 
Since U a (k) e SX(2,R), one obtains under the con- 
straints (34), (39) and (40) 



lim \U a (k)] 



12 



ad — 1 

7 



/3, 



(43) 



which can be shown in a direct manner by noticing 
Eqs.(30)-(32) and (35). Since v = 0(l/a 2 ), v + = 
0(1/ a) and w_ = 0(1 /a) respectively, a close exam- 
ination of the higher order terms in each equation is 
required to observe the cancellation among the terms 
of 0(l/a 2 ), 0(l/a) and O(a ). 

For the e potential of strength u, we have 



u 

72 : 



v (a) 
v + (a) = V-(a) 



2 

u 



(44) 
(45) 



which is exactly what we have shown in [18]. 

Up to now, we have excluded the case of 7 = 0. 
For this case, we start by assuming v to be constant 
(independent of the distance a), which is appropriately 
determined later. With this assumption, we see from 
Eq.(21) 



7 = a Hm o [W„(fe)] 21 



0. 



Moreover, Eqs.(22) and (23) are reduced to 

[U a (k)} n = 1+ (2a + 0(a 2 ))v+ + 0(a), 
[U a (k)} 22 = 1+ (2a + 0(a 2 ))v_ + 0(a), 

respectively. Eqs.(47) and (48) indicate that if 
a - 1 

v+(a) 



v-(a) = 



2a ' 
5-1 

2a 



(46) 



(47) 
(48) 



(49) 
(50) 



for small a, we attain 

a ]im [^( fc )]n = 
a }™ +0 [Ua(k)} 22 = S = 1/a. 



(51) 
(52) 



Here ad — 1 is imposed since det U a (k) — 1 and 7 = 0. 
It is easy to see from Eqs.(49) and (50) that the condi- 
tion aS — 1 is equivalent to 

v + + v- + 2v + v-a = 0. (53) 

Using Eqs.(33), (49), (50) and (53), we have 

, + + ,_ + [ZY a ( fe) ] 21 , + ,_ = ^- 1 4 ^- 1) 

+O(o). (54) 

We thus obtain the estimate of Eq.(24) for small a; 

[U a (k)] 12 = v + + v- + v + v (v + + w_)a 
+ [U a (k)} 21 v+v- +0(a) 

= a + 6 A +2 v + O(a). (55) 



Eq.(55) shows that if the constant vq is taken as 
4/3 



v = 



a - 



5 + 2 



(56) 



with an arbitrary real (i, one obtains the limit 

lim[U a (k)] 12 =p. (57) 

a >+0 

This completes the case of 7 = 0. In a particular case 
of a = 5 = 1 and j3 = v, Eqs.(49), (50) and (56) 
are reduced to = »- =0 and vq = v, namely a 
single 5 potential of strength v as expected. [Strictly 
speaking, Eq.(56) excludes the case of a = 5 = — 1; 
U = —Vs(—f3). However, this is easily realized by using 
a single 5 (v = —{3, v + = v_ = 0) together with the 
phase replacement from 9 to 9 + n.] 

Our findings in this section are summarized as fol- 
lows. The Hamiltonian H a in Eq.(17) produces all pos- 
sible connection conditions in Eq.(l) in the small a limit 
together with suitably renormalized strengths; In case 
of 7 4= 0, Eqs.(39), (40), (34) and (28) for v+, V-, v 
and A respectively, whereas in the case 7 = 0, Eqs.(49), 
(50) and (56) for v + , v_ and v together with Eq.(28) 
for A. 

3. Numerical Examples 

In order to show the validity of the three 5's approx- 
imation, we perform numerical tests in this section. 
For simplicity, we consider the bounded region [X\, X2] 
(Xi < < X2) together with the Dirichlet boundary 
condition <p(X\) = f(X2) = 0. The eigenvalues k n for 
the generalized point interaction are determined by 



[g(0,k n ;X 2 )Vg(0,k n ;-X 1 )} 21 =0 



(58) 
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Fig. 2 The solid line shows the wave function of the tenth (up- 
per) and the thirteenth (lower) eigenstate for the generalized 
point interaction with a = 3, (3 = —2, 7 = —7 and <5 = 5. The 
wave number is fcio = 0.894964 and k\z = 1.130869, respec- 
tively. For comparison, the approximated wave function with 
three <5's potential with the distance a = 0.2 is exhibited by the 
broken line. In this case, the wave number is fcio = 0.905264 
and fci3 = 1.142775, respectively. 

with Eqs.(2) and (12). The associated (not necessarily 
normalized) eigenfunction is calculated by 



<Pn(x) 



[G(0,k n ;x - Xi)} 21 , x<0, . . 

[e(O,ft n ;aOV0(O ) fc n ;--X-i)] 2 i. x > 0. l ° yj 



The eigenvalues for the three <5's approximation are de- 
termined by 

[G(0, K- X 2 - a)V a (A, k n )Q(0, k n ; -a - X ± )} 21 = (60) 

with Eq.(19), where A = 8/{2a) according to Eq.(28). 
In Eq.(19), the strengths of three S's are determined by 
Eqs.(34), (39), (40) for 7 4= and Eqs.(49), (50), (56) 
for 7 = respectively. The associated wave function is 
calculated in a similar manner to in Eq.(59). In the fol- 
lowing, we set the vector potential A = (8 = 0) since 
its effect is trivial. In the numerical calculations, we 
take the endpoints X2 = —X\ = 15.0 and the distance 
a = 0.2. 



(p 7 (x) 




Fig. 3 The solid line shows the wave function of the seventh 
eigenstate [kj = 0.775671 ) for the generalized point interac- 
tion with a = 5, (3 = 3, 7 = and 8 = 0.2. For comparison, 
the approximated wave function {kj = 0.775312 ) by three <5's 
potential with the distance a = 0.2 is exhibited by the broken 
line. 



Figure 2 shows the eigenfunction of two low-energy 
(tenth in upper and thirteenth in lower) eigenstates for 
the case of a = 3, (3 = —2, 7 = —7 and 6 = 5. The solid 
line is a calculation for the point interaction, whereas 
the broken line is the wave function approximated by 
three nearby S's. Figure 3 shows a typical example for 
the case of 7 = 0. We take a = 5, /? = 3, 5 = 0.2. In 
all cases, one can observe the wave function discontinu- 
ity around the point interaction, which is successfully 
reproduced in the corresponding eigenfunction approx- 
imated by three S's potential with appropriately renor- 
malized strengths. However, it is worthy to mention 
that the approximated ip n always has n — 1 nodes on 
the intervals as long as the distance a is kept finite, 
whereas this is not necessarily the case in the limit of 
a — ► +0, namely for the case of the generalized point 
interaction. Indeed, the lower part in Fig. 2 shows that 
ipi3 for the exactly point case loses two nodes in the 
vicinity of the interaction and as a result it has only 
ten nodes on the interval (X\, X2). On the other hand, 
the approximated (continuous) wave function behaves 
around the three S's in a somewhat complicated manner 
to ensure the relation between the number of nodes and 
the quantum number n mentioned above. Clearly, the 
convergence of the wave function in the small a limit 
is not uniform on the interval including the origin in 
general. 

4. Conclusion 

We have constructed one-dimensional generalized point 
interactions characterized by four parameters in the 
small distance limit of equally spaced three S's poten- 
tials. A constant vector field added on the interval be- 
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tween the two side (5's has no effect except for changing 
the phase of wave function. The remaining part, the 
transfer matrix U a {k) can be made converge into an 
arbitrarily fixed special linear matrix by adjusting the 
strengths of the three <5's according to the distance in 
an appropriate manner. Though a divergent term is 
inevitable for each strength in the small distance limit, 
it takes a remarkably simple and experimentally realiz- 
able form. Numerical examples support our model and 
show a satisfactory coincidence with the corresponding 
zero-range case even at a finite distance level. 
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